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IRI  PROJECT  REPORT  #4  -  FINAL  REPORT 

Reference  Number  F61 775-0 1-WE033 

Period  from  February  2003  until  August  2003  and  overview  over  the  entire  project 

1.  Summary  of  proposed  research 

1. 1  Description  of  services,  according  to  the  contract 

The  contractor  has  investigated  the  nonlinear  microwave  response  of  high-temperature 
superconducting  (HTS)  films  for  high-performance  microwave  components  that  are  used  in 
communication,  radar,  and  other  systems.  Contractor  has  provided  and  characterized  HTS 
films,  and  delivered  them  to  Dr.  John  Derov,  AFRL/SNHA.  Contractor’s  proposal  “Nonlinear 
microwave  response  of  HTS  microwave  devices”,  is  incorporated  as  attachment. 

1.2  Summary  of  research  objectives ,  according  to  the  proposal 

The  goals  of  the  project  were  a  comprehensive  understanding  and  control  of  the  non¬ 
linear  microwave  response  of  HTS  films,  in  terms  of  intermodulation  distortion  and  power- 
dependent  surface  impedance.  These  goals  have  been  achieved  to  a  major  extent,  even  though 
some  puzzles  related  to  the  microscopic  origin  of  the  nonlinearities  remain  to  be  solved. 

While  numerous  research  groups  and  high-tech  companies  have  demonstrated  HTS  de¬ 
vices  and  subsystems  of  unprecedented  performance,  the  nonlinear  response  of  the  supercon¬ 
ducting  components  eventually  limits  wider  application,  especially  in  narrowband,  sharp  cut¬ 
off,  receive  filters,  or  in  high-power  transmit  filters.  Among  the  basic  features  of  nonlinear 
microwave  response  are  a  RF-current-dependent  microwave  surface  impedance,  Zs(/r f),  and 
third-order  two-tone  frequency  intermodulation  distortion  (IMD).  Both  effects  have  been 
studied  for  superconductors,  with  emphasis  on  HTS  [1].  We  have  found  in  the  course  of  this 
project  that  the  dielectric  substrates  can  also  lead  to  undesirable  nonlinearities,  namely  a  RF 
field-dependent  loss  tangent,  tan SfEf),  and  IMD,  up  to  temperatures  close  to  Tc  [2,3]. 

The  microscopic  origins  of  IMD  in  YBaCuO  and  MgO  are  poorly  understood  at  present. 
Also,  little  is  known  about  the  relations  between  the  power  dissipation  over  microwave  cy¬ 
cles,  which  is  implicit  in  the  concept  of  a  surface  impedance,  and  the  IMD  caused  by  two 
closely  spaced  frequency  tones.  Investigating  and  modelling  the  nonlinear  behaviour  hence 
challenges  the  development  of  ultra-high  performance  front-ends  and  the  optimisation  of  new 
functional  microwave  materials  [3-6]. 
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1. 3  Development  of  the  project  and  topical  focus  of  the  final  half-year 

During  the  course  of  the  project,  in  total,  16  epitaxial  YBCO  films  on  MgO  were  pre¬ 
pared  in  5  batches  at  QinetiQ  Malvern.  Additional  MgO-substrates,  provided  by  STI,  LG 
Electronics  (Seoul,  Korea),  Lincoln  Laboratory,  and  eSCeTe  (Cr-doped  MgO),  were  used  for 
Nb  deposition  and  comparative  studies.  The  microwave  measurements  were  performed  at  LL 
and  AFRL  using  the  stripline  resonator  technique.  It  turned  out  that  measurements  on  unpat¬ 
terned  films,  which  were  intended  to  be  established  at  19  GHz  [4],  were  not  essential  for  the 
progress  and  success  of  the  project.  Therefore,  our  studies  focused  on  the  extended  experi¬ 
mental  body  available  from  LL  and  AFRL.  The  results,  analyzed  and  modeled  by  the  contrac¬ 
tor,  have  been  published  in  eight  papers  (partly  submitted  to  reviewed  journals  and  partly 
contributed  to  international  conferences)  [2,3,  7-12].  Our  results  and  observations  attracted 
attention  from  the  international  scientific  community  and  stimulated  research  in  different 
places.  The  project  has  helped  to  identify  routes  to  improve  the  power  handling  capability  and 
reduce  the  IMD  of  YBCO  films.  Further  potential  benefits  could  result  from  the  development 
of  novel  devices  based  on  the  nonlinear  dielectric  loss  tangent  of  MgO. 

During  the  last  phase  of  the  project,  we  have  studied  the  nonlinear  microwave  response 
of  a  batch  of  overdoped  YBCO  films,  part  of  which  were  additionally  subject  to  Ca-doping 
([7-9]  and  Sec.  3.1).  The  Ca-doping  was  performed  by  the  group  of  Professor  Jochen  Mann- 
hart  at  the  University  of  Augsburg  [13],  whose  contribution  we  gratefully  acknowledge.  The 
aim  of  this  study  was  to  find  out  how  far  the  improvement  of  the  linear  microwave  perform¬ 
ance  reaches  into  the  overdoped  range. 

Other  subjects  investigated  during  the  final  half-year  are  the  development  of  a  phe¬ 
nomenological  model,  which  relates  the  IMD  to  the  nonlinear  features  of  the  kinetic  conduc¬ 
tivity  of  the  superconductor. 

2.  Summary  of  achievements  during  the  final  half-year  period 

We  have  extended  the  correlations  between  the  DC  critical  current  density  of  supercon¬ 
ducting  YBCO  films  and  their  nonlinear  microwave  response  deeper  into  the  region  of  over¬ 
doped  films.  The  RF  critical  current,  and  hence  the  1-dB  compression  point,  as  well  as  the 
third-order  intercept  point,  continue  to  move  towards  higher  power  for  films  with  higher 
oxygen  content  and,  possibly,  ordering  of  the  oxygen  vacancies.  We  could  not  observe  any 
indication  of  saturation,  still  leaving  potential  for  further  improvement.  The  carrier  density 
can  be  adjusted  by  low-temperature  anneals  in  activated  oxygen,  or  by  Ca-doping  [13,14]. 
Oxygen  content  and  ordering  seem  to  be  the  key  parameters  to  optimize  the  linear  microwave 
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response  of  YBCO.  The  ordering  of  the  oxygen  vacancies  is  expected  to  be  affected  by  ad¬ 
justing  the  growth  stoichiometry  of  the  YBCO  films,  and  by  the  annealing  schedule  during 
and  after  film  deposition.  However,  reliable  procedures  remain  to  be  identified. 

We  have  developed  a  phenomenological  model  to  relate  the  IMD  signal  to  the  nonlinear 
contributions  of  the  surface  impedance  of  superconductors.  The  model  is  based  on  a  lumped- 
element  equivalent  circuit,  which  has  the  advantage  of  intuitive  access.  The  resulting  equation 
for  the  intermodulation  power  has  been  confirmed  by  a  separate  rigorous  analysis  of  inter¬ 
modulation  in  a  distributed  resonator  [15],  similar  to  the  derivations  in  [16,17].  By  comparing 
the  IMD  signal  predicted  by  the  model  with  our  experimental  data,  we  find  the  resistive 
nonlinearity  dominating  the  inductive  contribution  to  the  IMD,  in  contrast  to  the  theoretical 
expectations  from  [16].  Our  model  indicates  also  how  data  obtained  for  different  samples  or 
with  different  resonator  geometries  can  be  compared. 

In  addition  to  the  superconducting  nonlinearities,  we  have  also  observed  that  the  dielec¬ 
tric  loss  tangent  of  MgO  becomes  RF  field-dependent,  due  to  resonant  absorption  by 
Mg(OH)2  impurity  states.  The  dielectric  response  causes  marked  additional  microwave  losses 
at  temperatures  below  about  50  K.  The  dielectric  losses  include  contributions  from  power- 
independent  relaxation,  leading  to  pronounced  loss  peaks  around  25  K  in  the  2-GHz  region. 
Significantly  enhanced  intermodulation  is  another  consequence  of  the  microwave  loss  mecha¬ 
nisms  in  MgO.  Because  of  the  high  quality  of  our  latest  batch  of  superconducting  films,  we 
were  able  to  observe  the  influence  of  the  dielectric  nonlinearities  up  to  temperatures  of  about 
70  K.  The  dielectric  nonlinearites  are  stronger  in  more  slowly  cooled  films.  On  the  other 
hand,  full  oxygenation  and  good  oxygen  ordering  require  some  time  spent  at  relatively  low 
annealing  temperatures.  The  optimization  of  YBCO  films  on  MgO  for  microwave  device 
applications,  therefore,  requires  proper  adjustment  of  the  annealing  schedules  during  the 
different  phases  of  film  preparation. 

Beside  the  relevance  of  our  results  for  HTS  microwave  applications,  our  findings  could 
affect  any  high-frequency  application  involving  MgO,  or  similar  dielectrics.  The  effects 
brought  about  by  resonant  absorption  could  also  be  used,  in  principle,  for  new  types  of  signal 
detection  or  modulation  in  cryogenic  dielectric  microwave  devices. 
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3.  Nonlinear  microwave  response  of  YBCO 
3. 1  Properties  of  the  latest  batch  # 3318 

The  results  achieved  during  the  preceding  phases  of  the  project  have  shown  that  the 
nonlinear  surface  impedance  (i.e.,  resistance  and  changes  of  the  reactance)  as  well  as  the 
intermodulation  signal  decrease  for  increasing  oxygen  content  [4-6],  This  observation  let  us  to 
prepare  another  batch  of  overdoped  films,  including  the  approach  of  Ca-doping,  to  increase 
the  number  of  free  charge  carriers  as  much  as  possible  [8,13].  Fig.  1  sketches  the  parameters 
varied  in  batch  #5.  Table  I  presents  a  complete  list  of  the  films  studied  in  this  project. 

All  four  YBCO  films  of  batch  #5  were  annealed  at  540  °C  in  atomic  oxygen  after  depo¬ 
sition  (similar  to  the  treatment  of  #4B).  This  step  aimed  at  getting  a  comparable  and  -  suppo¬ 
sedly  -  homogenized  distribution  of  oxygen  vacancies.  The  40-nm  thick  cap  layers  on  #5A- 
#5C  were  then  deposited  at  the  University  of  Augsburg  by  pulsed  laser  deposition  at  tempera¬ 
tures  around  750  °C.  The  cap  layer  was  stoichiometric  YBCO  for  #5A  and  30%-Ca-doped 
YBCO  for  #5B  and  #5C.  The  cap  layers  of  #5A  and  #5C  were  removed  by  ion  milling  after 
deposition.  All  films  underwent  a  final  anneal  at  540  °C  in  atomic  oxygen  for  30  minutes. 

According  to  the  data  in  table  I,  all  films  with  cap  layers  displayed  Tc-values  well  below 
90  K,  with  #5B  and  #5C  having  the  lowest  Tc  <  86  K,  while  that  of  the  reference  film  #5D 
remained  close  to  90  K.  The  additional  deposition  and  corresponding  annealing  steps  reduced 
the  Tc-values,  probably  by  hole  doping  as  expected.  Interestingly,  though,  the  DC  critical 
current  densities,  which  still  exceeded  most  previous  JC6o-values,  decreased,  too. 

Fig.  2  compares  the  temperature-dependent  effective  surface  resistance  of  the  films  of 
batch  #5  with  a  number  of  films  from  other  batches.  The  most  striking  feature  is  the  pro¬ 
nounced  maximum  of  /?eff{T)  for  batches  #4  and  #5.  We  attribute  this  peak  to  dielectric  re¬ 
laxation-type  losses  in  the  MgO  substrate,  as  the  slow  cooling  after  the  deposition  of  the 
YBCO  base  film  could  have  favored  the  nucleation  of  Mg(OH)2  precipitates  [6].  Below  about 
20  K,  the  dielectric  losses  of  MgO  become  dominated  by  resonant  effects,  i.e.,  they  become 
power-dependent  and  hence  nonlinear.  We  expect  that  the  minimum  of  Retr  near  50  K  results 
from  the  power  dissipation  in  the  superconducting  film,  which  increases  as  Tc  is  approached. 
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Figure  1.  Schematic  description  of  the  films  of  batch  #5  [8]. 
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Batch/sample 

TC(K) 

Jc60 

(MA/cm2) 

c-axis  (nm) 

oxygen,  level 

remarks 

#  1  (2746) 

i  a 

89.7! 

12.96 

Slightly  underd. 

Damaged 

IB 

2.67 

1.1715 

Near  optimum 

As-grown 

1C 

2.23 

1.1716 

Near  optimum 

Molecular  02 

ID 

— 

1.1703 

Near  optimum 

Contact  anneal 

#2(2809) 

2A 

91.74 

11.40 

1.1703 

Optimal 

Qualitat.  data 

2B 

90.68 

5.82 

1.1725 

Atomic  oxygen 

2C 

84.18 

1.1686 

Overdoped 

Atomic  oxygen 

2D 

— 

— 

Damaged 

#3(2997) 

3  A 

90.54 

1.1734 

Near  optimum 

As- grown 

3B 

92.44 

6.88 

1.1724 

Optimal 

Damaged 

3C 

69.3 

2.67 

1.1752 

Underdoped 

Atomic  oxygen 

3D 

91.8 

1.1702 

Overdoped 

Atomic  oxygen 

#  4 (3321) 

4A 

87.75 

8.38 

n.m. 

Overdoped 

As-grown 

4B 

86.21 

11.00 

n.m. 

Overdoped 

Atomic  oxygen 

540  °C,  single 

4C 

86.75 

9.90 

n.m. 

Overdoped 

Atomic  oxygen 

400  °C,  single 

4D 

92.21 

5.60 

n.m. 

Overdoped 

Atomic  oxygen 

40(fC,  iterative 

#5(3318) 

5A 

86.79 

n.m. 

Overdoped 

YBCO  cap  removed 

5B 

82.84 

n.m. 

Overdoped 

Cacap 

5C 

85.59 

n.m. 

Overdoped 

Ca  cap  removed 

|  5D 

89.40 

14.60 

n.m. 

Overdoped 

without  cap  layer 

Table  I.  Updated  summary  of  characteristic  parameters  of  all  YBCO  films  on  MgO  prepared  at  QinetiQ  Mal¬ 
vern.  All  batches  were  nucleated  at  680  °C,  Y-rich,  and  grown  at  700  °C,  either  stoichiometric  (2746  and  2997), 
or  slightly  Cu-rich  (2809,  3318,  3321).  The  Tc-  and  Jc-values  were  derived  from  DC  magnetization  measure¬ 
ments  (further  remarks  see  text).  The  oxygenation  level  was  concluded  from  comparison  with  other  samples.  The 
cap  layers  used  for  batch  5  were  40  nm  thick  and  consisted  either  of  stoichiometric  or  30%-Ca-doped  YBCO. 

Although  the  microwave  losses  in  MgO  and  in  YBCO  have  different  temperature  scales,  it  is 
illustrative  to  compare  Ren(T)  for  batch  #5  on  a  T-scale  normalized  to  Tc  (Fig.  3).  Among 
batch  #5,  the  reference  film  #5D  displays  the  lowest  surface  resistance.  The  three  films,  which 
have  (had)  cap  layers,  display  enhanced  surface  resistance  at  intermediate  temperatures. 
Height  and  width  of  the  dielectric  loss  peak  vary  among  all  four  films,  indicating  that  the 
nucleation  of  defect  dipoles  in  MgO  is  strongly  temperature-dependent  in  the  range  of  the 
applied  deposition  temperatures,  i.  e.,  between  690  and  750  °C.  The  film  with  the  Ca-rich  cap 
on  it  (#5B)  displays  the  highest  surface  resistance,  which  could  result  from  deteriorated  mi¬ 
crowave  properties  in  the  region  of  high  Ca-concentration  [8]. 


Final  IRI  report  SPC01-403  3  3/02-8/03,  page  6  of  19 


Figure  2.  Temperature  dependence  of  the  total  microwave  losses  for  YBCO  films  on  MgO  from  four  different 
batches.  Batches  #1  and  #2  were  cooled  after  film  deposition  at  a  cooling  rate  of  10  °C/min,  batches  #4  and  #5 
at  the  much  slower  cooling  rate  of  5  °C/min.  The  slow  cooling  enhances  the  dielectric  losses  of  MgO. 


The  data  in  table  I  and  Fig.  3  contain  no  indication  of  the  beneficial  effects  from  Ca- 
doping,  which  rather  were  observed  in  granular  YBCO  films  [13,18].  Vice  versa,  the  #5-films 
are  concluded  to  have  very  few  large-angle  grain  boundaries,  if  any,  and  oxygenation  levels 
well  in  the  overdoped  region.  Both  conclusions  are  confirmed  by  Fig.  4,  which  displays  the 
/RF-dependence  of  the  surface  resistance  and  changes  of  the  surface  reactance  for  sample  #5  A 
in  comparison  with  the  three  films  of  batch  #3,  for  T  =  50  K  (table  I  and  [5]).  The  range  of 
linear  microwave  response  for  #5A  is  almost  twice  as  large  as  that  for  the  overdoped  film  #3D 


Figure  3.  Total  microwave  losses  for  the  films  of  batch  #5  versus  reduced  temperature. 
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Figure  4.  RF-current- dependence  of  the  surface  resistance  (solid  curves,  filled  symbols)  and  changes  of  the 
surface  reactance  (dashed  curves,  open  symbols)  for  four  films  covering  the  range  of  strongly  underdoped  (#3C) 
to  overdoped  (#5  A).  The  current-dependence  of  the  surface  resistance  is  stronger  than  that  of  the  surface  reac¬ 
tance  in  all  cases.  Both  types  of  nonlinearity  decrease  with  increasing  oxygenation  level. 

(see  also  comparative  data  in  Sec.  3.2),  yielding  a  critical  RF  current  of  7Cjrf  ~  2  A.  This  value 
corresponds  to  a  magnetic  flux  density  of  5rf  ~  60  mT,  which  comes  close  to  the  lower  criti¬ 
cal  field  expected  for  bulk  single-crystalline  YBCO  for  in-plane  shielding  currents.  The  linear 
range  is  expected  to  be  limited  at  this  field  level  by  vortex  nucleation  and  motion  [1].  For  the 
given  resonator  geometry,  with  a  resonant  volume  of  ~  0.003  cm'3,  the  power  circulating  in 
the  resonator  at  7c,rf  ~  2  A  amounts  to  7>cjrc  ~  170  W.  If  this  value  were  compared  with  a  disk 
resonator  at  the  same  frequency,  which  occupies  a  volume  larger  by  about  a  factor  of  210,  the 
scaled  value  of  PaK  would  be  36  kW.  Finally,  accounting  for  an  approximately  fivefold  en¬ 
hancement  of  the  RF  current  at  the  edges  of  the  stripline  compared  to  the  total  current,  the 
circulating  power  would  amount  to  ~  1  MW.  This  numerical  example  demonstrates  the  high 
film  quality  achieved  with  batch  #5,  which  comes  close  to  or  even  exceeds  the  Pcjrc-values 
achieved  with  TMoio  disk  resonators  at  this  temperature  [1,19]. 

We  noted  earlier  [5],  that  the  temperature  dependence  of  the  resistive  nonlinearity  can 
be  scaled  by  the  RF  current,  independent  of  the  oxygenation  level.  We  expect  a  similar  sca¬ 
ling  behaviour  for  the  surface  reactance.  The  critical  RF  current,  therefore,  is  a  major  parame¬ 
ter  describing  the  nonlinear  surface  impedance.  Altogether,  we  observe  a  systematic  and 
reproducible  variation  of  7c>rf,  at  a  given  temperature,  by  a  factor  of  ~  7  for  the  investigated 
range  of  oxygenation  levels  from  strongly  underdoped  to  highly  overdoped  (c.f..  Fig.  7). 

Fig.  4  also  shows  that  the  resistive  nonlinearity  is  stronger  than  the  reactive,  i.e.,  indue- 
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output  power  at  fundamental  (dBm) 


Figure  5.  Evolution  of  the  IMD-signal  with  temperature  for  sample  #5A.  The  enhanced  IMD  due  to  the  nonlinear 
dielectric  loss  tangent  of  MgO  disappears  at  elevated  temperature.  The  power  dependence  of  this  evolution  is 

depicted  in  greater  detail  in  Fig.  13. 

tive,  nonlinearity  in  our  high-quality  YBCO  films  on  MgO.  This  is  an  important  experimental 
fact,  as  numerous  groups  tend  to  neglect  the  resistive  contributions  to  nonlinear  effects  like 
intermodulation  distortion  or  harmonic  generation,  especially  in  resonant  devices  (Sec.  3.3). 

Fig.  5  illustrates  the  temperature  dependent  IMD  signal  versus  output  power  at  the  fun¬ 
damental  tones  for  #5A.  The  broad  plateau  below  20  K  is  associated  with  the  dielectric 
nonlinearity  in  MgO  [4].  For  increasing  temperature,  the  IMD  signal  decreases  continuously 
at  low  power,  but  passes  through  a  minimum  at  intermediate  power  levels.  Cuts  through  the 
diagram  at  constant  power  are  shown  in  Fig.  13  below  and  discussed  in  further  detail  there. 

Fig.  6  compares  the  IMD  levels  of  the  films  introduced  in  Fig.  4.  The  data  clearly  show 
decreasing  IMD  for  increasing  oxygenation  level.  At  intermediate  output  power,  we  find  a 
variation  of  the  IMD  level  among  the  four  samples  of  ~  30  dB.  According  to  general  argu¬ 
ments  and  confirmed  by  our  model  in  Sec.  3.3,  a  third-order  nonlinearity  like  intermodulation 
distortion  depends  quadratically  on  the  quadratic  coefficients  of  Rs(Irf)  and  AYs(/rf).  Since 
the  nonlinear  surface  impedance  can  be  scaled  by  the  current,  the  IMD  signal  is  expected  to 
depend  on  the  fourth  power  of  1//c,rf-  An  increase  of  /Cjrf  by  a  factor  of  7  corresponds,  there¬ 
fore,  to  a  reduction  of  the  IMD  level  of  ~  34  dB,  which  is  close  to  what  we  observe.  This 
consistency  indicates  that  the  two  types  of  nonlinearity,  Zs(/rf)  and  IMD,  are  caused  by  a 
single  mechanism,  despite  the  very  different  time  scales  involved.  We  are  also  led  to  specu¬ 
late  that  the  dependence  of  7c,rf  on  the  oxygenation  level  bears  a  clue  to  the  microscopic 
understanding  of  the  nonlinear  microwave  response  of  YBCO. 
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output  power  at  fundamental  (dBm) 

Figure  6.  Comparison  of  the  IMD  signal  versus  output  power  at  the  fundamental  tone  for  the  four  films  described 
in  Fig.  4,  for  T=  50  K.  The  intermodulation  distortion  decreases  for  increasing  oxygenation  level. 

3.2  Updated  correlation  plots 

The  trends  for  7Cjrf  and  Pimd  discussed  along  with  Figs.  4  and  6  are  generalized  in  Figs. 
7  and  9,  to  include  all  samples  studied  until  now.  Fig.  7  displays  the  correlation  between  the 
RF  critical  current  and  the  DC  critical  current  density  for  T=  50  K.  The  relationship  is  almost 
linear  with  no  sign  of  saturation.  The  highest  critical  currents  were  achieved  with  the  latest 
batch,  #5.  The  graph  yields  a  mean  slope  5/c>rf/3Jc  ~  10'7  cm2.  Interpreting  this  value  as  an 


Figure  7.  Correlation  between  the  RF  critical  current  and  the  DC  critical  current  density  for  T~  50  K.  /c  RF 
increases  almost  linearly  with  Jc,  with  a  slope  ~  0.1  A/(MA/cm2).  The  four  films  of  batch  #5  (table  I)  display  the 
highest  critical  parameters  of  all  samples  studied  until  now. 
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effective  cross-section  of  the  current  path,  whose  height  is  approximately  limited  by  the  film 
thickness,  t  ~  350  nm,  we  find  a  range  of  widths  between  20  and  75  pm.  If  the  critical  current 
densities  at  DC  and  RF  frequencies  were  equivalent,  this  length  scale  would  indicate  a  typical 
domain  size  of  undisturbed  current  flow.  It  is  worth  noting  that  this  dimension  appears  consis¬ 
tent  with  the  distribution  of  the  RF  current  across  the  width  of  the  stripline  (150  pm).  RF  mea¬ 
surements  on  different  resonator  geometries  could  be  used  to  verify  the  equivalence  of  RF  and 
DC  critical  currents. 

According  to  Fig.  7,  the  dominant  source  of  the  nonlinear  microwave  response  has  the 
same  origin  as  the  breakdown  of  the  inductive  shielding  of  low-frequency  currents.  It  has 
been  argued  that  twin  domain  boundaries  or  oxygen  vacancies  could  reduce  the  critical  cur¬ 
rent  density  by  depinning  [20],  This  point-of-view  would  be  consistent  with  the  correlation 
between  the  DC  critical  current  density  and  the  c-axis  parameter  of  epitaxial  YBCO,  shown  in 
Fig.  8  for  some  samples.  However,  the  scatter  in  this  and  other  correlation  plots  indicates  that 
the  electronic  properties  of  the  YBCO  films  are  not  determined  by  the  oxygen  content  alone, 
but  likely  by  the  ordering  of  the  oxygen  vacancies,  too.  This  aspect,  which  could  be  investi¬ 
gated,  e.g.,  by  systematic  studies  of  the  orthorhombicity  of  the  films,  should  be  a  major  focus 
of  research  activities  following  this  project. 

Fig.  9  displays  the  linear  relation  between  the  third-order  intercept  (TOI),  relative  to  the 
output  power  at  the  fundamental  tones,  and  the  DC  critical  current  density.  The  overall  varia¬ 
tion  of  the  TOI  is  slightly  less  than  the  30  dB  expected  from  the  IMD  data  in  Sec.  3.1.  This 
discrepancy  can  be  attributed  to  deviations  of  the  IMD  slope  m  from  the  exact  value  3,  as 


lattice  parameter  c-1 .16  nm  (pm) 

Figure  8.  Correlation  between  the  DC  critical  current  density  at  T-  50  K  and  the  c-axis  lattice  parameter  for 

some  films  studied  in  this  project  (c.f.,  table  I). 
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Figure  9.  Correlation  between  the  third-order  intercept  at  T=  20  K  and  the  DC  critical  current  density  at  T=  50 
K.  The  TOI  increases  almost  linearly  with  Jc,  with  a  slope  —  2  dBm/(MA/cm2).  The  film  #1B  falls  apart  from  the 

overall  trend  for  reasons  not  yet  known. 


illustrated  by  Fig.  10.  The  slopes  observed  for  the  latest  batch  of  films,  m  ~  2.8,  are  higher 
than  could  be  anticipated  on  the  basis  of  the  preceding  results  [5].  The  majority  of  films  now 
displays  slopes  between  2.5  and  3.  In  retrospective,  it  seems  as  if  the  small  IMD  slopes  of 
#2C  and  the  films  of  batch  #4  could  have  been  due  to  some  kind  of  defects.  At  least  for  batch 
4,  this  assumption  is  in  accordance  with  the  steep  inductive  nonlinearity  at  low  RF  currents 
[6],  which  pointed  to  the  presence  of  weak  links  that  may  have  been  the  result  of  damage 
introduced  during  the  patterning  process. 


2 

critical  current  density  Jc(50K)  (MA/cm  ) 


Figure  10.  Correlation  between  the  IMD  slope  at  T- 40  K  and  the  DC  critical  current  density  at  7"=  50  K. 
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As  for  the  nonlinear  surface  impedance,  we  conclude  that  the  intermodulation  distor¬ 
tion,  too,  is  determined  to  a  major  extent  by  the  critical  parameter  7c,rf,  or  Jc,  respectively. 

3.3  IMD  model 

We  have  developed  a  theoretical  model  to  provide  a  solid  base  for  the  analysis  of  our 
IMD  measurements  described  above.  While  other  groups  have  treated  similar  problems  in  the 
literature  [16,17,21],  our  aim  was  to  relate  the  measurable  IMD  signal  to  other,  measurable, 
quantities  like  the  quadratic  coefficients  of  the  current-dependent  resonant  frequency  and  un¬ 
loaded  quality  factor.  As  a  starting  point,  we  have  used  the  equivalent  circuit  depicted  in  Fig. 
11a.  This  lumped-element  representation  was  successfully  employed  for  the  calculation  of  the 
linear  current  in  a  resonator  for  a  given  input  power  [22],  We  have  adapted  this  model  to  the 
nonlinear  case  as  indicated  by  Fig.  lib.  There,  Fimd  is  the  voltage  generated  by  the  nonlinear 
inductor  and  resistor.  We  have  used  the  following  ansatz  for  the  nonlinear  circuit  elements: 

X(Irf)  =  X0+X2I2rf,  (1) 

where  X  stands  for  the  reactance  or  resistance  per  unit  length,  a>L  or  R.  The  quadratic  coeffi- 


Figure  1  la.  Equivalent-circuit  model  of  a  capacitively  coupled  resonator  (red)  terminated  by  a  voltage  source  at 
the  input  and  a  matched  load  at  the  output.  This  circuit  is  used  for  the  calculation  of  the  linear  current  [22]. 


Figure  1  lb.  Equivalent-circuit  model  used  for  the  calculation  of  the  IMD  current  [15].  This  circuit  is  based  on 
the  all-parallel  version  of  the  linear  circuit  depicted  in  Fig.  1  la.  The  voltage  source  contains  contributions  to  the 

IMD  current  from  both  inductor  and  resistor. 
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cient  X2  is  in  units  of  Q/A2m;  I  is  the  total  current  at  the  two  input  frequencies  (fundamental 
tones): 

/  =  70  •  (cos  co,/  +  cos  co20 .  (2) 


We  have  checked  the  results  of  our  calculation  with  a  separate  analysis  of  a  distributed  reso¬ 
nator  [15],  similar  to  the  treatment  of  Dahm  and  Scalapino  (DS)  [16].  The  analysis  of  lumped- 
element  and  distributed  resonators  led  to  the  same  expression  for  the  IMD  signal  Pimd,  to 
within  a  factor  rc/3,  which  we  neglect  in  the  following: 


PIMD(dBm)  =  -91 .74  +  40  •  log[rv(l  -  rv)Q0]+  3  Pin(dBm)  +  10- 


(3a) 

(3b) 


1/2 

with  Zo  =  47  Q  the  characteristic  impedance  of  the  stripline  resonator,  rv  =  (P0ut/Pm)  the 
voltage  insertion  ratio,  go  the  unloaded  quality  factor,  Pin  and  Pout  the  powers  at  input  and 
output  of  the  resonator,  and  fi  and  Q2  the  quadratic  coefficients  of  resonant  frequency  and 
unloaded  quality  factor.  The  latter  are  related  to  J/2  by: 


[2  _  Pi  Qi _ P^_ 

f0  24’  g0~  v 


The  contribution  of  the  inductive  nonlinearity  to  go  =  coLIR  is  smaller  by  a  factor  1/go  than  its 
contribution  to  /2  and  can  hence  be  neglected. 

Eq.  (3)  contains  the  expected  cubic  dependence  of  Pimd  on  Pm,  which  is  an  inherent 
feature  of  the  model.  Pimd  depends  also  on  the  g-factor  of  the  resonator,  via  go  and  rv.  If  the 
quality  factor  itself  depends  on  the  input  power,  the  functional  dependence  PiMD(Pin)  will  be 
modified,  and  the  original  behaviour  of 4  and  Q2  can  be  masked.  The  result  in  Eq.  (3)  resem¬ 
bles  the  DS-result  except  for  a  different  prefactor  (3/2)  that  essentially  results  from  a  different 
normalization,  and  the  different  weighting  of  resistive  and  inductive  nonlinearities.  We  have 
convinced  ourselves  that  both  types  of  nonlinear  contributions  are  equally  important.  Ways  to 
access  their  relative  weight  by  experiment  have  been  proposed  in  [15,17]. 

In  order  to  account  for  the  explicit  g-dependence  of  Pimd  in  Eq.  (3)  properly,  the  equa¬ 
tion  can  be  rewritten  in  terms  of  the  loaded  g-factor  and  the  g-factor  associated  with  the 
(symmetric)  coupling,  grad: 


(5b) 
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Eq.  (5a)  relates  the  IMD  signal  to  the  output  power  at  the  fundamental  tone,  whereas  Eq.  (5b) 
refers  to  the  circulating  power,  ECirc  =  PoJQkA.  It  is  worth  noting  that  the  Q-scaling  of  Pimd 
involves  Q\2  in  both  cases. 

We  have  adapted  our  model  to  calculate  the  IMD  signal  expected  from  fi  and  Qi  for  the 
two  films  #3D  and  #5A,  with  the  results  summarized  in  Fig.  12.  The  quadratic  coefficients 
were  extracted  from  the  corresponding  RF-current  dependences  by  fitting  polynomials  of 
order  2  (quadratic)  or  4  (quartic)  to  the  data.  We  find  satisfactory  agreement  between  calcula¬ 
tion  and  measurement  at  elevated  power.  We  attribute  the  discrepancy  at  low  power  to  a 
weakness  related  to  the  determination  of  the  quadratic  coefficients,  as  these  cannot  be  sepa¬ 
rated  properly  from  the  linear  background. 

For  both  samples  studied  in  Fig.  12,  the  resistive  nonlinearity  was  stronger  than  the  in¬ 
ductive  contribution  to  the  IMD  signal,  as  summarized  in  table  II.  This  table  also  illustrates 
how  sensitive  the  extracted  nonlinear  coefficients  depend  on  the  way  how  they  are  deter¬ 
mined.  We  note  that  the  relative  strength  of  the  resistive  nonlinearity  is  smaller  for  the  better 
film  #5A.  Further  studies  should  investigate  the  weighting  of  the  two  nonlinear  contributions 
in  greater  detail,  e.g.,  by  phase-sensitive  measurements  [23]. 


|  Type  of  fit  f 

#3D 

|  #5  A 

[  quadratic  | 

6.82 

1  3.18 

|  quartic  f 

33.33 

j  16.45 

Table  II.  Summary  of  results  obtained  for  the  ratio  of  the  normalized  quadratic  coefficients,  (Qi/2Qo2)/(f2/fo), 
deduced  from  quadratic  and  quartic  fits,  respectively,  as  explained  in  the  text. 


Figure  12.  Comparison  between  the  measured  IMD  signal  (symbols)  and  the  calculated  results  (lines)  for  sam¬ 
ples  #3D  (blue)  and  #5A  (red)  for  T=  50  K  using  quadratic  (solid  lines)  and  quartic  fits  (dashed  lines).  The 
agreement  at  elevated  power  results  from  a  dominant  contribution  of  the  resistive  nonlinearity  (table  II). 
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4.  Nonlinear  microwave  response  of  MgQ 

The  consequences  of  resonant  absorption  in  MgO  have  been  analyzed  and  discussed  in 
detail  in  previous  reports  [4-6].  We  focus  here  on  new  results,  which  could  be  identified 
thanks  of  the  high  film  quality  of  batch  #5.  Fig.  13  displays  the  temperature  dependent  mi¬ 
crowave  losses  for  #5A  for  two  extreme  power  levels  (dark-blue  and  light-blue  curves,  left- 
hand  ordinate).  The  two  most  striking  features  are  the  minimum  of  the  losses  around  T=  40  K 
at  a  relatively  high  level,  Ret r  ~  10  pQ,  and  the  splitting  into  low-power  and  high-power 
branches  below  T  =  20  K,  which  is  associated  with  the  resonant  absorption.  We  attribute  the 
minimum  to  the  superposition  of  two  counter-acting  mechanisms:  the  decrease  of  dielectric 
losses  caused  by  defect  dipole  relaxation,  and  the  increase  of  losses  in  the  superconductor  as 
the  critical  temperature  is  approached.  We  have  found  that  the  losses  due  to  dielectric  relaxa¬ 
tion  were  higher  in  YBaCuO  films  that  were  more  slowly  cooled  from  their  deposition 
temperature  to  ambient  (c.f.,  Fig.  2  and  Fig.  14).  This  is  consistent  with  the  nonlinear 
behaviour  of  the  substrate  being  due  to  resonant  absorption  by  Mg(OH)2  impurity  states  [24]. 

Fig.  13  also  shows  temperature  dependent  IMD-data  for  constant  output  power  at  the 
fundamental  tone  (dark-green  and  light-green  solid  curves,  right-hand  ordinate).  These  data 
represent  "vertical  cuts"  through  the  curves  of  Fig.  5,  for  two  selected  values  P0ut  =  -60  dBm 
and  -40  dBm.  Both  curves  feature  a  pronounced  drop  of  IMD  at  the  temperature  where  the 
resonant  absorption  in  MgO  disappears.  Interestingly,  for  the  lower  output  power,  the  IMD- 
signal  keeps  dropping  up  to  the  highest  temperature,  T~  70  K,  while  the  IMD-level  increases 


H 


Figure  13.  Temperature  dependence  of  Reff  (left-hand  ordinate)  at  the  lowest  measured  power  (light-blue)  and  at 
its  minimal  value  (dark-blue).  Displayed  on  the  right-hand  ordinate  is  the  IMD  signal  for  two  levels  of  output 
power  (light-green:  Pout  =  -60;  dark-green:  -40  dBm).  The  dotted  curves  represent  the  IMD  signal  after  scaling 
with  respect  to  the  temperature  and  power-dependent  loaded  Q-factor  [Eq.  (5)]. 
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smoothly  above  about  50  K  for  the  higher  output  power.  This  increase  is  due  to  the  growing 
influence  of  the  superconducting  nonlinearity  near  Tc.  However,  at  70  K,  the  IMD  at  low  and 
high  power  still  only  differ  by  about  four  orders  of  magnitude,  not  the  60  dB  that  would  be 
expected  from  a  purely  cubic  power  dependence. 

These  features  of  the  temperature-dependent  IMD  could  be  affected,  at  least  partly,  by 
the  power-  and  temperature-dependent  (7-factor  [see  Eq.  (5)].  To  check  for  such  ambiguities, 
we  have  scaled  the  measured  IMD  signal  by  [0l(2OK)/£?l(T)]2,  leading  to  the  green  dotted 
curves  in  Fig.  13.  The  comparison  of  solid  and  dotted  curves  reveals  that  the  loaded  quality 
factor  does  affect  the  measurable  IMD  signal,  and  that  g-scaling  is  necessary.  But  it  shows 
also  that  the  main  features  discussed  above  remain  unchanged.  This  result  confirms  the  en¬ 
hanced  IMD  caused  by  the  nonlinear  loss  tangent  of  MgO. 

In  our  earlier  data,  the  nonlinearities  caused  by  the  superconducting  films  masked  the 
dielectric  nonlinear  contributions  from  the  substrate  already  at  temperatures  below  about  50 
K.  Our  new  results,  which  were  obtained  for  films  with  a  very  weak  superconducting  nonlin¬ 
earity,  reveal  that  the  dielectric  losses  in  MgO  persist  up  to  T-  70  K.  This  result  needs  to  be 
taken  into  account  for  most  microwave  applications  of  HTS,  which  are  tailored  to  a  tempera¬ 
ture  region  between  50  and  80  K.  Unfortunately,  using  our  superconducting  stripline  tech¬ 
nique,  we  were  not  able  to  measure  the  nonlinear  loss  tangent  and  IMD  of  the  MgO  alone, 
which  could  persist  to  even  higher  temperatures.  We  suggest  dielectric  resonators  to  study  the 
nonlinear  response  of  MgO  in  more  detail.  We  note  in  passing  that  the  dielectric  nonlinearity, 


E0(V/m) 

Figure  14.  Parametric  plot  monitoring  the  dielectric  loss  tangent  and  the  saturation  field  of  MgO  for  T~  4  K.  The 
plot  contains  data  for  YBCO  films  (red  colour)  as  well  as  for  Nb  reference  films  (blue).  The  effect  of  slow 
cooling  in  terms  of  enhanced  ARttt  and  tendentially  increased  E0  is  clearly  visible.  The  dashed  lines  are  guides  to 

the  eye  (c.f.,  [6]). 
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which  is  undesirable  for  microwave  applications  of  HTS,  may  be  beneficial  for  RF  signal 
processing  with  dielectric  materials. 

Fig.  14  summarizes  the  parameters  describing  the  nonlinear  loss  tangent  of  MgO,  for  all 
samples  studied  in  the  course  of  this  project.  The  results  for  batch  #5  added  new  evidence  that 
slow  cooling  of  the  films  after  deposition  are  detrimental  to  the  dissipation  losses  in  MgO. 
According  to  the  data  for  #4D,  the  enhanced  loss  tangent  can  be  reduced,  at  least  partly,  by 
post-deposition  annealing.  This  possibility  opens  the  potential  for  simultaneous  optimization 
of  the  oxygenation  state  of  the  superconducting  films  and  the  loss  tangent  of  the  MgO. 

5.  Open  questions,  perspectives  for  future  research  activities 

The  aims  of  the  project  specified  in  Sec.  1  could  be  achieved  to  a  great  extent.  However, 
there  remain  interesting  questions  to  be  solved,  which  could  be  investigated  in  future  research 
activities  following  this  project. 

5. 1  Nonlinearities  in  the  superconductor 

1 .  Study  a  batch  of  films  rapidly  cooled  after  film  deposition,  and  with  the  oxygenation  state 
optimized  by  post-deposition  annealing. 

2.  Study  homogeneously  Ca-doped  films  for  a  set  of  Ca-concentrations. 

3.  Investigate  the  orthorhombicity  of  the  YBCO  films  used  for  RF  measurements. 

4.  Perform  RF  measurements  using  stripline  resonators  of  different  geometry  to  investigate 
the  correlation  between  /Cjrf  and  Jc  in  greater  detail. 

5.  Perform  phase-sensitive  IMD  measurements  to  judge  the  relative  weighting  of  resistive  and 
inductive  nonlinearities. 

6.  Investigate  the  causal  relations  underlying  the  microscopic  mechanism  of  the  nonline¬ 
arities:  superconducting  order  parameter  -  oxygenation  state  -  Jc  -  nonlinear  microwave 
response. 

7.  Refine  studies  of  third-harmonic  generation  (THG).  Our  IMD  model  (Sec.  4)  can  be  used  to 
analyze  THG-effects,  too.  Because  the  third-harmonic  is  non-resonant  in  the  stripline  resona¬ 
tor  due  to  mode  dispersion  [10],  the  magnitude  of  the  measured  THG  signal  depends  strongly 
on  the  frequency  offset  between  the  third  overtone  and  the  third  harmonic.  Since  studies  of 
the  THG-signal  could  help  distinguish  between  resistive  and  inductive  nonlinearities,  special 
stripline  resonators  could  be  envisaged  where  both  frequencies  of  interest  are  resonant. 


Final  IRI  report  SPC01-4033  3/02-8/03,  page  18  of  19 


5.2  MgO-anomalies 

8.  Study  the  influence  of  annealing  and  quenching,  to  identify  optimal  conditions  for  a  low 
loss-tangent. 

9.  Perform  measurements  of  the  microwave  response  of  bare  MgO  by,  e.g.,  using  a  dielectric 
resonator). 
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